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Abstract 

In choosing a family of histories for a system, it is often convenient to choose 
a succession of locations in phase space, rather than configuration space, for 
comparison to classical histories. Although there are no good projections onto 
phase space, several approximate projections have been used in the past; three 
of these are examined in this paper. Expressions are derived for the probabil- 
ities of histories containing arbitrary numbers of projections into phase space, 
and the conditions for the decoherence of these histories are studied. 
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I. INTRODUCTION 



A great deal of work has been done recently on the use of the decoherence formalism 
to describe quantum mechanical systems |T|-^|. These systems can be described in terms of 
decoherent histories, which can be assigned probabilities obeying classical probability sum 
rules. While, in principle, a history could be described in terms of any set of variables, 
most of the work has focussed on histories of particles in configuration space. The simplest 
and most fine-grained such history is just the classical trajectory of a particle, specifying its 
exact position at every moment in time. Such histories do not decohere, however. Instead, 
one must consider considerably coarse-grained histories, in which a position is given only 
at certain discrete times, and only within certain finite intervals. A history can then be 
specified by a string of indices a^, stating which interval the particle is in at time t{. 

Another important class of histories, though, would be descriptions of a system as being 
in cells of phase space at successive points in time. A small amount of work has been done 
on this subject P-13|, but they have not been tackled in full generality. 



A. The Decoherence Functional 

The decoherence functional is a functional on pairs of histories of a quantum mechanical 
system. One simple description of the functional (though not the most general) has the 
form: 

D[a,a'} = Tr{p a " n (t n ) • • ■ P^tJpP^U) • ••/%(*„)}. (1.1) 

In this expression, p is the initial density matrix of the system. The -P^(ti) are Heisenberg 
projection operators onto Hilbert space. At each time ti these projection operators represent 
different alternatives for the system. In terms of Schrodinger projections P, these time 
dependent projections can be written P(t) = e lHt / h p e - tHt / h . A complete set of histories has 
an exhaustive set of alternatives at each time, 
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E^(**) = i- 



(1.2) 



A particular choice of the {ctj} represents one particular history, which we will denote a for 
brevity. Thus, saying that a given history a occurs implies that alternative cti occurs at 
time ti, a 2 at time t 2 , and so forth. 

The probability of a given history occuring is just given by the diagonal elements of 



In order for these histories to obey the classical probability sum rules, we must require that 
the set of histories decoheres. The usual requirement for this is that the off-diagonal terms 
of the decoherence function vanish, 



This is actually a sufficient, but not a necessary condition for decoherence. All that is truly 
required is that the real parts of the off-diagonal terms vanish. Most physically decoherent 
systems, however, display this stronger form of decoherence; we will see this in the cases 
that we consider. 



The most common type of problem treated at present is that in which the variables 
are divided into a system and a reservoir, or environment. In this case one traces over the 
reservoir variables, and is left with a reduced density matrix on only the system variables. 
Instead of the simple time evolution operator e ~ lHt ^ h , the system evolves according to a 
somewhat more complicated transition matrix or propagator T. In terms of path integrals 
this is 



D[a,a'}: 





D[a, a'] = 0, a ^ a'. 



(1.4) 



B. The Transition Matrix 



T(xf, x tf] Xi, x, 



■;,«=/ 
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where the integral is over all paths x(t) and x'(t) from ij to tf which begin at Xj and x\ and 
end at x / and x'j, respectively, ^[^(t)] is the action of the system variables independent 
of the reservoir, and W[x(t),x'(t)} is the Feynman- Vernon influence phase arising due to 
the interactions with the reservoir ||14|| . The model most commonly considered is one that 



was developed in the study of Brownian motion [HJ, in which a one-dimensional particle 
described by a single variable x interacts with an infinite bath of harmonic oscillators via a 
linear or weakly nonlinear potential, but this formalism is quite general. 
The reduced density matrix p(x; x') evolves straightforwardly: 

^.j^M^,,^^^. (,6) 

Thus, the decoherence functional can now be written 

D[a,a'} = Tr.{i£T(. • ■ T(P>P a \) • • (1.7) 

With the projections P being onto intervals of coordinate space, it is very easy to write 
the decoherence functional as a constrained path integral over x and x' . For phase-space 
projections, the form of the decoherence function is more complicated, as we shall see. While 
there are no true projectors onto cells of phase space (as there are for intervals of coordinate 
space), there are a number of approximate projectors, and we shall consider these one at a 
time. 

In dealing with phase space, it is natural to consider other representations of the density 
matrix, most obviously the Wigner distribution: 

w(X,p) = - f dee^V(X + e/2;X-e/2). (1.8) 

7T J 

This distribution acts in many ways like a probability distribution in phase space, with the 
major exception that it can be negative in localized regions. The time evolution of w(X,p) 
is also described by a transition matrix: 

w(Xf,p f )\t f = J J dXidpi T w (Xf,pf,t f ;Xi,pi,ti)w(Xi,pi)\ ti . (1.9) 



T w {X f ,p h t f -X h p h U)^ J J dtfdti e^f-MM* 

xT(X f + f - + |,X, - |,t,), (1.10) 

where this T is the same transition matrix defined above ( |1.5|) . We shall see that the 
expressions for the probabilities of phase space histories are described very naturally in 
terms of Wigner distributions. 



II. PROBABILITIES OF PHASE SPACE HISTORIES 



As mentioned above, there are no true projections onto cells of phase space ||10|| . This 
is essentially a side-effect of the uncertainty principle, which prevents both x and p from 
being localized simultaneously. However, for cells larger in area that h, we can approximate 
projections reasonably well. 

In making these calculations, we will find that it is useful to work in terms of w(X, P) 
and T w . For this we use the inverses of (|1.8| ) and ( |1.10| ), namely 

p{x;x') = J dp e- ip{x - x ' )/n w{{x + x')/2,p) (2.1) 

and 
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xT w ((x f + x' f )/2,p f ,t r ,(x t + x'J/2 }Pi} U). (2.2) 



A. Consecutive X and P Projections 

While there are no good projections onto phase space cells, projections onto intervals 
in coordinate or momentum space are perfectly well-defined and straightforward. They are 
just 

i-Xi+Ax/2 

P x = I dx \x)(x\, (2.3a) 

Jxi-Ax/2 
m+Ap/2 

Pn= \ , dp\p)(p\. (2.3b) 



We can imagine using a projection to determine which interval of x the system is in, followed 
a short time later by a projection onto an interval of p. If we let the time between these two 
"measurements" go to zero, we can make use of the relation 

(x\p) = e ipx/h . (2.4) 

This is unsatisfactory in a number of ways. The most obvious is that x and p are non- 
commuting variables, so that P x and P p are also non-commuting. The order in which one 
makes these measurements matters, particularly if the intervals are fairly small (compared 
to h). If we are interested in fairly large cells in phase space, this is of less importance; for 
classical and quasiclassical systems this is often the case. 



Measurements of this type were treated by Halliwell [10]. He also considered another 
type of two-projection measurement: a pair of successive position measurements, separated 
by a small time interval At, with the momentum determined by the time of flight between 
x\ and x-i- I have not considered this type of measurement, as it is ill-defined as At — > 0, 
and therefore requires non-trivial consideration of the system's time evolution between the 
two position projections. For a system with complex dynamics this is difficult. 

The exact projections ( ^.3[ ) used above are less convenient for the purposes of calculation, 
though they are certainly more correct. For ease of computation, therefore, it is customary 
to use approximate Gaussian projections, 

P- x = —I— r dx e-( x -* i)2/Ax2 \x)(x\, (2.5a) 



y^Ax J-oo 

P*= j~Jp e-^ 2 ^ 2 \p)(p\. (2.5b) 

Using these projections, the probabilities just reduce to a product of Gaussian integrals, 
which can, with a little algebra, be easily solved. Using the expressions ( [2.1| ) and ( |2.2| ) for a 
history with iV measurements of cells of phase space centered on (x\,pi), (22,^2), • • • (xn,Pn), 
we get 
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xW(X , p )5( Zl - (X + ^/2))5{z[ - (X - Co/2)) 

( Zl - x,) 2 + (z[ - x,f (h - Pl ) 2 + (k[ - Pl ) 2 



x exp 



Ax 2 



Ap 2 



+jr[ik , 1 z[ - k lZl ) + X 1 (k 1 - k[) + ^i(pi - fci/2 - fci/2) - £ Po 
xT(X 2 ,j9 2 ,t2;Xi,pi,ti) 



(2.6) 



x 



The integrals over z, z', k, k', and £ are all simple, and yield 

Pa= ( 2 j)N-l J d{ X i}d{Pi}™( X 0,Po) 



x exp 



x exp 



2 2 Ax 2 

(Pi ~ Pi) 2 ~ — ( X « " ^i) 2 - ^ 

Ax 2 



Ap 2 ^ 1 ^' Ax 2 
xT(X 2 ,p 2 ,t 2 ;X 1 ,p 1 ,t 1 ) 



Ap 2 
2h 2 



T ( P o - Pi) 2 - ^(Xo - X,) 2 



2 (ps - Pi) 2 ^ -w iP2 ~ nf ~W {X2 ~ Xs)2 



Ap 2 ^° Ax 2 
xT(X 4 ,p 4 ,t 3 ;X 3 ,p 3 ,t 2 ) 



x 



x exp 



{X 



Ax 2y " 2N ' 2 2[h 2 + Ax 2 Ap 2 /4]. 

Note that the expression for the probability behaves very reasonably, i.e., the evolution 
after a "measurement" continues to be centered about the measured values of X and p, with 
a spread determined by the size of the phase space cell. 

A measurement of p followed by a measurement of X produces an expression very similar 
to the above, and is readily evaluated by the same methods. The differences are chiefly 
notable when the cell size is small compared to H. 



x N y 



Ax 2 (p 2N -2 ~Pn) 



(2.7) 



B. Coherent State Projections 

The closest thing to a true projection onto a cell in phase space is probably the coherent 
state projection \p,x)(p,x\ centered on (p,x). While these are true projections, they are not 
orthogonal: 
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\(p, x\p', x')\ 2 = exp[-(x - x'f/2a - a(p - p'f /2% 



(2.8) 



Also, these states are overcomplete. Thus, phase space histories built from coherent states 
cannot be truly decoherent, and can only be even approximately decoherent if a discrete 
sample of them (e.g., the states corresponding to a lattice of points in phase space) is taken. 
In a coordinate basis we can represent a coherent state as 

1 



(x\p, x) 



(trt) 



1/4 



exp 



(x — x) 2 



2a 



+ 



ipx 



(2.9) 



This expression is useful in evaluating the probability of a coherent state history. 

When we consider a history of N "measurements" in phase space using coherent state 
projections we get ciii expression £lii£l1cl gous to ( |2.7D , which can (again) be solved for the 
probability: 

4^ 



2na 



d{Xi}d{pi} w(X ,p ) 

(Xq-xx) 2 (Xx-xi) 2 {po-pifo- (Pi-Pi) 2 ^ 



x exp 

L a a 

xT(X 2 ,p 2 ,t 2 ;X 1 ,p 1 ,t 1 ) 



h 2 



K 2 



XT \X 2 N-2, P2N-2, £jV) -^2^-3,^2^-3) ^JV-l) 
(X 2 Ar_2 - X N ) 2 (p 2N -2 ~ Pn) 2 CT 



x exp 



a 



n 2 



(2.10) 



The general behavior of the probabilities is very similar to that in the first case we 
considered, but even cleaner and easier to see. Coherent states are an excellent way of 
representing phase space histories. 

There is one other kind of approximate projection that we could consider. It is not, in 
my opinion, a very attractive one, but it has been used in the literature, and so might as 
well be treated here. Consider approximate "projections" of the form 



— - — / dpdx \p, x)i 
ipAx J 



-(p-p) 2 /Ap 2 -(x-x) 2 /Ax 5 



(p,x\. 



ixApA 

If we consider iV measurements of this form, the probability becomes 



(2.11) 



S 



Pa = 2a + **A? + oAp 2 Ax 2 ' d W d <*> 



x exp 



(X - xi) 2 + (Xx - x x ) 2 (po - pO 2 + (px - p^ 2 



o'- 



er + Ax 2 h 2 + oAp 2 

2 if- 2 A 2 , _ A „2 A „2 /*2\ / (^0 ~ X X ) 2 (p - p x ) 



-(o 2 Ap 2 /h 2 + Ax 2 + oAp Ax /% ) v u ; x/ 9 , + , 
v y 1 vi ! \2o(o + Ax 2 ) 2(h 2 + oAp 2 ) 

xT(X 2 ,p2,t 2 ;X 1 ,p 1 ,t 1 ) 



xT(X 2 N-2, P2N-2, tN] X 2 Ar_3,p 2 JV-3? *JV-l) 

0{0 + Ax 2 ){p 2N -2 - pN? + + oAp 2 ){X2 N -2 ~ X N f" 



x exp 



(2.12) 



2oh 2 + a 2 Ap 2 + Ax 2 ft 2 + aAp 2 Ax 2 
Again, the same sort of qualitative behavior, but a much uglier expression. 

III. DECOHERENCE OF PHASE SPACE HISTORIES 



While the above expressions are highly intuitive in their qualitative behavior, we have (in 
a sense) been putting the cart before the horse. It is meaningless to assign a probability to a 
history without first being assured that the set of histories described is decoherent. There is 
nothing in the expressions above to prevent one from choosing extremely tiny cells in phase 
space, with areas small compared to %\ yet such histories are certainly not decoherent, as 
they flagrantly violate the uncertainty principle. 

Unfortunately, while we can write expressions for the probabilities without having to 
know much about the physics of the system (i.e., the actual behavior of the transition 
matrix T), in order to actually calculate them, or to say much about decoherence, we need 
to know something about the path integrals. 

Except in the case of quadratic systems, these integrals are not exactly solvable. Limited 



treatments of this case have been considered elsewhere [|T^,[n|. Most interesting systems, 
however, include nonlinearities. This can be handled in one of three ways: numerically; in 
perturbation theory; or in the semiclassical limit, where solutions are peaked about the "clas- 
sical trajectory." The first approach is robust, but does not lend itself to general arguments. 
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While almost any approach will eventually have to be treated numerically to calculate ac- 
tual values for probabilities or decoherence functional elements, one would hope to get a 
rough idea as to a systems behavior before invoking that numerical machinery. The second 
approach, perturbation theory, is the most commonly adopted. When the nonlinearities are 
weak, the path integrals can be approximated with considerable precision. Unfortunately, 
many interesting cases (e.g., chaotic systems) cannot be treated in this fashion; for them, 
their nonlinearities are intrinsically important. The last approach is limited to systems with 
sufficient mass and inertia to resist quantum fluctuations [[§. This is useful in considering 
either the classical limit of quantum systems, or in estimating quantum effects in otherwise 
classical systems, and is the approach we will adopt here. 

Earlier work has concentrated on distinguished systems interacting with a large reservoir 
or environment whose degrees of freedom can be neglected. As has been shown, such systems 
give rise to decoherence functionals with probabilities peaked about classical trajectories. 
The transition matrix for such a system has the form (in the limit of a large thermal reservoir) 

T(Xi, £1, ti\ X Q , £ , to) = 

J 5X5£ exp (MX(t) + dV/dX(X(t)) + 2MjX(t) - g(t)^£(t)dt 

+ ^ir^ C em + M ^ _ + 0(e3) }' (3,1) 

where X and £ are variables defined by 

X=^(x + x'), (3.2a) 
i=x-x', (3.2b) 

and the reservoir temperature is T. This is basically a toy system, consisting of a single 
one-dimensional particle of mass M moving in an arbitrary potential V(X). The interaction 
with the reservoir provides the dissipative term and a thermal noise; it is this noise which 
causes the system to decohere. |£| is a measure of how far "off-diagonal" the decoherence 
functional is; for large |£| it will clearly be strongly suppressed. This is "medium-strength" 
decoherence as defined by Gell-Mann and Hartle f2j. 
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Since large £ is suppressed, we can neglect the higher order terms in £ with good accuracy. 
This makes the £ path integral purely quadratic, and therefore solvable. Doing this integral 
yields 

T(Xi, £1, ti] X , £ , *o) — 

+(i/h)M(£ 1 X 1 -SoXo)}, (3.3) 

which is clearly peaked about the solution to the classical equation of motion 

MX(t) + dV/dX(X(t)) + 2M 1 X{t) = g(t), (3.4) 

more and more strongly in the limit of large M. 

Let X d (t) be the solution to the above classical equation with the boundary conditions 
X c \(to) = X and X c \(ti) = X±. We can then define a new variable rj(t) 

V (t)=X(t)-X d (t). (3.5) 

Clearly r)(t) has boundary conditions i](ti) = i](t ) =0. As M becomes large, we can treat 
r)(t) as a small deviation, and approximate the path integral as 

T(X 1 , £i, ti] X , £ , to) = 



6V -p{-^£ 1 (M^(t)+rfV/rfX 2 (X cl (t))^(t) + 2M 7 ,)(t)) 2 ^ 

+(i/h)M(^r )l - ^fj ) + ( i /h)M^ 1 X cl (t 1 ) - toX cl (t ))}. (3.6) 

This path integral is quadratic in rj and therefore solvable, at least in principle. This 
principle runs into a few problems in practice. It assumes that you know X c \(t) as a function 
of the boundary conditions. This is true only in very simple cases. In chaotic cases, it may be 
difficult to determine this function even numerically. Also, this integral contains (in essence) 
4th derivatives of rj, which complicate the calculation in some ways. Still, by making a few 
assumptions about the behavior of X c \(t), we can still extract some useful information from 
this expression. 
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Since for the purposes of determining decoherence we are really only interested in the £ 
dependence of T, it is straightforward, albeit tedious, to show that 

T{Xi, £1, ti, X , £ , to) = 
^exp|-^^(A 1 £ 2 + A 2 £o6 + A 2 £ 2 )+(z/fr)M(£ 1 X cl (t 1 ) - £ ^ci(*o))}- (3.7) 

Note that Aj = Xi(Xi, X , t\, to) and K = K(Xi, X , ti, t ). These functions are not es- 
pecially easy to calculate, but can be computed numerically if necessary. Simple calcula- 
tions along those lines seem to show that \i/{t\ — to) is relatively constant for (t\ — to) 
short compared to the dynamical time of the system and long compared to the decoher- 
ence time, at least for high-probability paths. For longer times, comparable to the dynamic 
timescale of the system in question, the Aj vary enormously in magnitude; numerical re- 
sults showed a variability of more than four orders of magnitude, though most results for 
Ai/(ti — t ) and A 3 /(ti — t ) clustered around certain values, and never became negligibly 
small. Ai£q + A 2 £o£l + ^3^i is, in any case, always a strictly non-negative quantity. For 
details of these calculations, see the Appendix. 

If we make the (admittedly highly questionable) assumption that the Aj are roughly 
constant for constant (t± — to), then we can estimate the level of decoherence achievable with 
phase-space projections. For simplicity, we will only look at the projections at a single time: 

Tr{.-.T(Pi(^T(---)P4(ti))--j (3.8) 



A. Consecutive X and P Projections 

Using the approximate X and P projections described in (|2.5b|) above, we can examine 



decoherence by looking at the off-diagonal elements, where the projections are centered 
on (x,p) and (x',p') respectively. A single pair of projections at time ti will multiply the 
decoherence functional by a factor 

dX 1 dX 2 d£ 1 d£ 2 dpdp'T(X 3 , £ 3 , t i+1 ; X 2 , f 2 , U) 

12 



, (X, - 6/2 - x'f {p-pf {p'-p'f 

X GXp s 



Ax 2 Ax 2 Ap 2 Ap 2 

+i(p - p')(X 2 - Xl ) + l -(p + p')(& - 6)} 

xT^.a.^Xo^o,^-!) (3.9) 



J dX 1 dX 2 K 2 esxp^-(^ 



(Xx - xf + (X x - x') 2 + (Ap 2 A:r 2 /2ft 2 )(X 2 - X 1 f 



Ax 2 

^_ ( (8h 2 M-fkT\ + h 2 Ap 2 ) (x - x') 2 



+ 



—i 



\2Ax 2 {8h 2 M-fkTX + h 2 Ap 2 + 16(M 7 A;TA) 2 Ax 2 + AM^kT\Ap 2 Ax 2 ) 

(h 2 + 2M<ykT\Ax 2 )[(p 2 - p) 2 + (p 2 - p') 2 } 

StfM^kTX + k 2 Ap 2 + !6(M<ykT\) 2 Ax 2 + AM 1 kT\Ap 2 Ax 2 

{2M 1 kT\Ax 2 ) [(pi - p) 2 + (pi - p') 2 ] 

8h 2 M<ykT\ + h 2 Ap 2 + l6(M<ykT\) 2 Ax 2 + AM 1 kT\Ap 2 Ax 2 

(ft 2 + AM^kTXAx 2 )^ -p') 2 - Ap 2 {p 2 - Pl ) 2 

2{8h 2 M 1 kT\ + ft 2 Ap 2 + 16(M 7 fcTA) 2 Ax 2 + AM 1 kT\Ap 2 Ax 2 ) 
(AHM^kTX) (x - x') (p + p' - 2 Pl ) 



8h 2 MjkT\ + h 2 Ap 2 + 16(M-fkT\) 2 Ax 2 + AM^kT\Ap 2 Ax 2 
i(p-p>){X 2 -X l 



h 



}, (3.10) 



where p\ and p 2 are MX d for boundary conditions {X(ti-i) = X , X(ti) — Xi} and {X(ti) = 
X 2 ,X(t i+1 ) = X 3 } respectively. 

A formidable expression indeed! One can, with difficulty, see that in general this factor 
will be suppressed for off-diagonal terms. If we simplify matters by taking the semiclassical, 
high-temperature limit, an examination of the real terms of the exponent show that for 
\x — x'\ 2 ~ 5x 2 this expression is suppressed by a minimum factor of 

5x 2 (\ h 2 



f teVl Tf \1 

6XP l A^ 2 \2 ~ AM^TXAx 2 ))' 



AM-fkTXAx 2 

A similar examination of the real p terms gives no similar comfort, for we find that there 
the minimum level of suppression is none at all! This doesn't mean that histories with 
differing p's do not decohere; the last imaginary term in the exponent oscillates extremely 
rapidly, and will tend to suppress all off-diagonal terms as X\ and X 2 are integrated over. 
This will work, in general if \p — p'\ 2 ~ 5p 2 is large compared to Ap 2 and h 2 / Ax 2 . To 
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suppress the off-diagonal terms of the decoherence functional by a factor e, where e 1, we 
must have 

5x 2 > -Arc 2 hie, 
5p 2 > -Ap 2 In e, 
Ax 2 Ap 2 > k 2 . (3.11) 



B. Coherent State Projections 

The results from coherent state projections are similar, but somewhat cleaner and easier 
to see. In this case the factor from the projections at one time U goes as 

\X 1 - xf + (X! - x'f + (X 2 - x) 2 + {X 2 - 



J dX 1 dX 2 K 2 exp^-(^- 



+ 



2a 

% 2 {x-x'f 



2h 2 a + STiM^kTXa 2 

(Pi ~ P? + (Pi ~ P' f + (P2 ~ P? + fa - gT - (P - pO 2 
2^ 2 /a + 8hM-fkT\ 

ih(x - x')(p 2 -pi) i(p -p')(X 2 - Xi) 



h 2 + AM^kTXa + h 



}• (3-12) 



Here again, we see that in the semiclassical limit this reduces to a minimal level of 
suppression 

5x 2 n n 2 



cxp 



a V2 SM-fkTXc 
for the x terms, and that bp 2 must be large compared to Ap 2 > h 2 /a. 



IV. CONCLUSIONS 



While no set of approximate phase space projections treated in this paper is completely 
satisfactory, they do serve to illustrate certain traits that phase space histories should pos- 
sess. Highly discontinuous trajectories are suppressed, and as one goes to the semiclassical 
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limit the probabilities of histories become peaked about the classical solutions. While pre- 
cise statements about decoherence are hard to make, given the difficulty of solving the 
problem for highly general systems, rough arguments indicate that the size of phase space 
cells needed for decoherence is much larger than that naively indicated by the uncertainty 
principle (AxAp ~ %). 
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APPENDIX: PATH INTEGRAL FOR THE TRANSITION MATRIX 

The path integral in ( |3.6|) is somewhat unusual in that it has only two boundary condi- 
tions (77(^0) — v(ti) — 0) for an integrand with four derivatives! Thus, the usual prescription 
for solving quadratic path integrals is not immediately applicable. 

This procedure can still be used, however, by the simple expedient of imposing two 
more boundary conditions, f](to) = Vo and rj{t\) = v%, and solving the path integral, then 
integrating the result over all values of vo and v±. 

The path integral to be solved is then 



ACKNOWLEDGMENTS 




F(ii,* ) exp{S d + i(M/h)(Z lVl - f t>o)}, 



(Al) 



where 
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(A2) 



e(t) = M(n + 2 7 >j + f(t)r)), (A3) 

m = < A4 > 

and F(ti,to) is an undetermined multiplier independent of the boundary conditions of rj. 
(Of course, since f(t) is defined in terms of X c i(t), this whole solution is dependent on X 
and X\. This dependence is complicated, as we will see.) 

The classical action S c \ is the action of the path rj(t) that obeys the classical equation of 
motion. For 



rti 

S = / L(r],rj,ri)dt 
J t 



the Euler-Lagrange equation is 



dt 2 \drj J dt \drj J \dr) J 
Plugging in our definitions for S[r](t)} and e(t), we get the equations 



i) + 2 7 ?) + f(t)r}= e{t)/M, (A6a) 
e - 2 7 e + f{t)e= 0. (A6b) 



The first of these follows from the definition ( |A3| ) , the second from 
The solution to the first equation is 

e(t) = Me^(A ai (t) + Ba 2 (t)), (A7) 

where a\ and a 2 are two independent solutions of the equation 

a = { 1 2 -f{t))a. (A8) 

This equation is not easily solved analytically in most cases. For slowly varying f(t) one 
can approximate the solution 
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a(t) « exp{± f ^ 7 2-/( s )ds} ; (A9) 

in any case, ( |A8| ) is readily solvable numerically. While any independent boundary conditions 
will work for a x and a 2 , a convenient choice is 

Ol (0) = l,Oi(0) = 0; 

a 2 (0) = 0,03(0) = 1. (A10) 
This is then plugged into the equation for 77 to give the solution 

v (t)=e^ t (Ca 1 (t) + Da 2 (t)) 



1 f a 1 (t)a 2 (s) - a 2 (t)a 1 (s) 



Jt \ai{s)a 2 {s) — a 2 (s)ai(s)/ 
Imposing the boundary conditions on 7] then gives us a set of equations involving A, B, C, 
and D: 

C = 0, (A12a) 

D = v , (A12b) 

u a 2 (*i) + ^(ai(ti)/ 2 - a 2 (t 1 )/ 1 )+ J B(a 1 (t 1 )/ 3 - a 2 {h)I 2 ) = 0, (Al2c) 

vod^h) + A(di(ti)/ 2 - d 2 {t l )I x )+B{d 1 {t 1 )h - d 2 {t x )I 2 ) = v 1 e rtl , (A12d) 

where 

j a = r ai ^) 2e27 ^ t (Ms) 

A, a'i(t)o 2 (t) — d 2 (t)ai(t)' 
/2= r«i a 1 (t)a 2 (t)e 2 ^dt 
it a'i(t)a 2 (t) - d 2 (t)ai(ty 

h= /" "ff^ (A15) 
/to d 1 (t)a 2 (t) - a 2 (*)ai(£) 

Again, these integrals are usually only solvable numerically. 
Solving for A and /? then yields 

^_ (di(ti)a 2 (ti) - ai(ti)d 2 (ti))-?3^o + e 7 * 1 (ai(ti) J 3 - a^Q/g)^! 
(ai(ti)d 2 (t!) - di(ti)o2(ti))(J 2 2 - A/3) 
= Aov + ^iWi, (A16) 
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B _ -(ai(ti)a 2 (ti) - ai(ti)d 2 (ti))/2^o - e 7 * 1 (ai(ti)/ 2 - «2(ti)-fi)^i 
~~ (ai(ti)d 2 (fi) -ai(ti)a 2 (t 1 ))(/ 2 2 - 7 1 7 3 ) ~~" 

= B v + B lVl . (A17) 

Plugging these results back into the definitions of e(t) and S[r](t)} gives a value for the path 
integral 



exp{s d +(iM//i)(£iUi - ^o)} 
M 



exp{-^(A 2 / 4 + 2ABI 5 + B 2 h) + i{M/%){^v x - ^o)}, (A18) 



where 



J 4 = [ 1 ai (t) 2 e 2 ^dt, (A19) 

/ 5 = [ 1 ai (t)a 2 (t)e 2 ^dt, (A20) 

J 6 = Z' 1 a 2 (t) 2 e 2 ^dt. (A21) 

7 1 



Clearly, the exponent is quadratic in t> and i>i. Integrating over these two boundary- 
conditions, we perform two gaussian integrals, and arrive at our final result 

J dv dv! J 5r) exp{%(*)]+i(M//i)(fiUi = 

Xexp{-^^(A 1 e 2 + A 2 66 + A 3 £ 2 )}. (A22) 



where 



_ A\h + + B\h 



A2 " 2(A 5 1 -A 1 5 ) 2 (/|-/ 4 / 6 ) ' (A23b) 



A 2 / 4 + 2A /j / 5 + B 2 h r . 9 _. 
A3 " 4(^-^0)^/1-/^)' (A23C) 



# = , ^ F(t 1; t ). (A23d) 

M^/(A 1 B -A B 1 y(hh-l!) 
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Clearly, these quantities are dependent on X and X\\ \ = \i(X , Xi, to, h) and K = 
K(X , Xi, t , ti). In principle, K can be calculated, but in practice it is not necessary to do 
so in order to make arguments about decoherence. It would be necessary to do so in order 
to actually compute the probability of a history. 

The derivations in this section, tedious as they are, can nevertheless be readily auto- 
mated. Once one determines the classical trajectory X c \ which corresponds to the boundary 
conditions X and Xi, determining the A« numerically is straightforward. I have used this 
technique to examine the values of the Aj for the forced, damped duffing oscillator model. 
For long times (ti — 1 ) the Aj varied enormously in magnitude as a function of X and X±. 
For times short compared to the dynamical time of the system, however, the Aj were nearly 
linear: \/(ti — t ) ss 0.085. For some boundary conditions they might become considerably 
larger, but in what was, admittedly, not an exhaustive sampling, none got much smaller. 
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